Anomalous Josephson current through a spin-orbit coupled quantum dot 
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For a general model of a mesoscopic multi- level quantum dot, we determine the necessary condi- 
tions for the existence of an anomalous Josephson current with spontaneously broken time-reversal 
symmetry. They correspond to a finite spin-orbit coupling, a suitably oriented Zeeman field, and 
the dot being a chiral conductor. We provide analytical expressions for the anomalous supercurrent 
covering a wide parameter regime. 

PACS numbers: 74.50.-|-r, 74.78.Na, 71.70.Ej 



Introduction. — The Josephson effect, where an equilib- 
rium supercurrent flows through a junction between two 
superconductors held at phase difference (f>, is of funda- 
mental importance in condensed matter physics, quan- 
tum information science, microelectronic applications, 
and metrology It has recently attracted renewed in- 
terest in mesoscopic and nanoscale junctions after the 
experimental demonstration of gate-tunable Josephson 
currents through junctions with just a few relevant elec- 
tronic levels ( "quantum dot" ) in a variety of material 
systems, e.g., InAs nanowires the 2D electron gas in 
semiconductors 0], and carbon nanotubes Q. One im- 
portant novel aspect arises because the spin-orbit interac- 
tion (SOI) strength a due to structural and bulk inversion 
asymmetries is often significant 0, 0], and theoretical 
work has therefore started to address SOI effects on the 
Josephson current in such devices 0, M, M, [13, [IHj [13, [I3l ■ 
So far this activity has mainly focused on 0- and tt- 
junctions (positive or negative critical current /c, respec- 
tively). A remarkable prediction concerns the possibility 
for an anomalous supercurrent la, flowing even at zero 
phase difference (</> = 0) if both a (suitably oriented) Zee- 
man field b and the SOI are present [i,[i3,[il. For </> = 0, 
the Hamiltonian is invariant under time-reversal symme- 
try (TRS) even when a 7^ and b 7^ 0, and the anoma- 
lous supercurrent thus spontaneously breaks TRS, which 
otherwise enforces /(</)) = — /(— <^) and hence /a = [l|. 
Anomalous supercurrents were first predicted in uncon- 
ventional superconductors [l^ . but were never observed 
there. In the tunneling limit, where the conventional 
current-phase relation (CPR) reads I{(f)) = Ic sin </>, this 
is equivalent to a, phase shift 0o, i.e., I{4>) = /csin(04-0o) 
and thus la = /cSin(/>o. The phase shift in such a "i^o- 
junction" could be observed in a SQUID containing one 
0- and one (/)o-junction via the shift of the diffraction 
pattern, or as spontaneous current in a superconducting 
ring containing a (/)o-junction. Both effects are tunable by 
external gate voltages (affecting the SOI), Zeeman fields, 
and by an orbital magnetic fiux. Moreover, a (^n-junction 
can also act as a superconducting rectifier [12l |. 

Mesoscopic systems contacted by conventional s-wave 
BCS superconductors could then yield a new class of 
systems with spontaneously broken TRS, and therefore 



exhibit anomalous supercurrents. Recent works have 
started to address this point. First, Ref. Q consid- 
ered a long ballistic one-dimensional Rashba quantum 
wire, where /a 7^ is tied to the Zeeman effect and to 
the difference between the velocities of right- and left- 
moving electrons. However, for physically realizable a, 
the reported la values turn out to be extremely small, 
/a (X a*, or are most likely inaccessible in experiments. In 
a mainly numerical study [T^ ]. the anomalous Josephson 
effect was also found for a multi-channel spin-polarizing 
quantum point contact. Finally, Buzdin [isj reported 
/q 7^ in junctions containing a noncentrosymmetric fer- 
romagnet as a weak link. There is clearly a need to 
systematically classify all ingredients necessary to ob- 
serve the anomalous Josephson effect. In this work, 
we compute the Josephson current through a generic 
phase-coherent mesoscopic system (an arbitrary multi- 
level quantum dot). Analytical predictions for la are 
provided for a physically important parameter regime. 
In addition, these results are supported by numerics. A 
necessary condition for la ^ emerges from our study: 
the quantum dot must be a chiral conductor, see Eq. 
below. This requirement was irnplicit in Ref. [3] but is 
apparently violated 15[ in Ref. [1J| . 

Model and exact solution. — As a generic model for the 
Josephson current through a quantum dot, we consider 
H = Hl + Hft + Ht + Hd, with two identical s-wave 
BCS superconductors {Hj^L/ji) of gap A held at phase 
difference and tunnel-coupled {Ht) to the quantum 
dot. With lead fermion operators Ci^/j^ k^ for spin a =t, j 
and momentum k, the BCS Hamiltonian reads (we often 
put h = I and, for simplicity, take the zero-temperature 
limit) 



For a = and b = 0, the closed dot can be described 
in terms of real- valued eigenfunctions Xn{f) [with r = 
{x,y,z)] for eigenenergy e„, where n = 1,...,N labels 
the relevant dot orbitals. Using fermion operators dna for 
these orbitals and disregarding strong correlation effects, 
the dot Hamiltonian in the presence of SOI and Zeeman 
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field can be written in the form (the PauH matrices (Tx,y,z 
act in spin space, spin indices are kept impUcit, and we 
absorb the magnetic factor gfJ-s/^ into b) 



H 



D 



E4 



b ■ cr] d„ 



cr d„ 



(1) 



We consider both the Rashba and Dresselhaus SOI due 
to structural and bulk inversion asymmetry, respectively, 
with the electric field pointing in z direction. The SOI is 
encoded in the antisymmetric matrices a^^y forming the 
vector [3 



ann' ^ — dr x„(r) 

TO ' 



sm{9)dx + cos{9)dy 



cos(6i)a^ - sin(e')ay ) 

(2) 

The overall SOI strength is represented by the character- 
istic inverse length a, while the angle 9 = Q corresponds 
to a pure Rashba and = 7r/2 to a pure Dresselhaus case 
0. The tunneling Hamiltonian is 



j — L I R,k7ia 



(3) 



where we make the inessential assumption of k- 
independent tunneling matrix elements. It is useful to 
define the Hermitian N x N hybridization matrices (in 
dot level space). 



- j — L I H.nn' 



(4) 



where v is the normal-state density of states in the leads. 
We find that the anomalous supercurrent is most pro- 
nounced if the Zeeman field points along a direction de- 
fined by the SO angle 9, namely h — B{cos9, — s'm9, 0)"^. 
When the magnetic field is orthogonal to this, the anoma- 
lous supercurrent is absent. The above model ignores or- 
bital magnetic fields. For the important case of a 2D dot, 
this is justified by choosing an in-plane magnetic field as 
above. 

Integrating out all fcrmionic degrees of freedom, the 
Josephson current for arbitrary system parameters is ob- 
tained: 



I{<P) = -— / duj trlnS 
h Jo 



(5) 



with a 4 A'' x AN matrix S{uj). Using auxiliary Pauli 
matrices t^,^, the matrix S can be expressed as [lOj 



S = + 



A 



X [{Tl + Tn) cos(0/2)ct,t, + {Tl - Tr) sin(0/2)aj,] , 

where E = diag(ei, . . . , cat) with e„ = e„ — o? jlm con- 
tains the dot level energies. The Zeeman field and the 
SOI are encoded in the w-independent AN x AN matrices 



with the vector of real anti-symmetric matrices [l6| 

^cos6'[l - 2cos(26') sin^(ax)] 
d^XudyXn' \ sin0[l -f 2cos(20) sin^(aa:)] 
cos(20) sin(2ax) 

(7) 

The matrices B^^y.z arc real and symmetric, and for the 
above Zeeman field given by 

Bnn' ^ B dr XnXn' 



(8) 



COS cos^(aa;) — cos(30) sin^(aa;) 
X I — sin0cos^(Q!a;) — sin(30) sin^(ax) 
— cos(26') sin(2ax) 



One can verify from the above expressions that no 
anomalous supercurrent can exist in the absence of ei- 
ther the SOI or the Zeeman field. 

Analytical approach. — We now derive the anomalous 
supercurrent la in the most relevant limit of weak SOI 
and weak Zeeman field, where la oc aB is small. More- 
over, the derivation below assumes that the off-diagonal 
entries in the F^/^ matrices Q are small against the 
diagonal entries. This condition is met in most cases 
of practical interest, e.g., if one dot level is resonantly 
coupled to the leads and all other levels are only weakly 
coupled, or when quasi-random phase shifts between dif- 
ferent tjn have to be taken into account. We then con- 
sider the limit — > 0, where under the above conditions, 
it makes sense to write S = Sq + Si in Eq. The 
"leading" part is diagonal in dot level space. 



Sn = -i 



A2 



(9) 



with To = diag(ri + r^). Similarly, using Ti = Tl 
Tr — To , the "perturbation" part is 



(0/2)A(rL - TR)<jy + [-iuj + AaxT,]Ti 
Vuj^ + A2 



(10) 

with Z given in Eq. ([6]). The anomalous supercurrent 
then follows by expansion of the tracelog in Eq. ([S]), 



2e 



E 



(-1)' 



(11) 



Using Eq. ^ , straightforward but lengthy algebra shows 
that both the n = 1 and n = 2 contributions always 
vanish. The leading contribution to la then comes from 
n = 3, where the part oc Fi in Si, see Eq. ([T0| . does not 
contribute at all. In the end, we arrive at the analytical 
expression 



^eA^ 



2 roo 



duj 



trj [TR,TL]-D-\b- A) 



Bx(Ty)Ty 



(6) 
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with the diagonal matrix D = uj'^ + E'^ + TI+ ^ Jl^^a Tp ■ 
The trace operation tr^ extends over dot level space only. 
Since b-A cx aB, we see that la oc aB as expected. Apart 
from having a finite SOI and an appropriately oriented 
Zeeman field, an additional condition must be satisfied 
in order to have 7^ 0: 



(13) 



which implies chirality for transport through the quan- 
tum dot. From numerical studies of the full current 
([5|), we find that Eq. is a necessary condition for 
-fa 7^ for arbitrary other parameters, and hence is 
not restricted by the conditions under which Eq. (fT^ 
has been obtained. As a consequence, a single-level dot 
(N = 1) can never allow for an anomalous supercurrent, 
and at least two relevant orbitals are required. Finally, 
we remark that Eq. (jl2p can be further simplified for 
max(r„„) ^ A, where we obtain 



la = ^trrf {{Tl 



(14) 

Basic explanation. — Having established that sponta- 
neously broken TRS is indeed possible in this system, see 
Eq. (|12p . we now give an intuitive physical argument to 
understand the origin of this effect. Consider the transfer 
of a Cooper pair through a quantum dot with N = 2 lev- 
els ioi (f> ~ 0, schematically shown in Fig.[TJ The Cooper 
pair has amplitude tL->_R {tn^i,) for transfer from left to 
right (right to left). For simplicity, we put = and as- 
sume real tunneling amplitudes ij„ in Eq. To lowest 
order in aB, SOI and Zeeman field combine to the term 
H' = {iAx+Bx)<Jx, with 2x2 matrices in dot level space, 

'^l ) ' ^^'^^'^ A cx a, and Bj. ~ ^ ^ q ^ 

To lowest order in the tj^ and in H' , processes like the 
one sketched in Fig.[lja) generate a contribution to tL_+_R. 
Specifically, the process in Fig. [TJa) yields 



StL^B, = {tLitm) {tLi iAB i_R2), 



(15) 



where the |-electron causes the first factor while the "f- 
electron switches levels via the product of A^ and Bx 
processes. In effect, the SOI and the Zeeman field con- 
spire to produce an effective orbital magnetic field via 
such processes, which then explains the broken TRS. 
With the group velocity w, the anomalous supercurrent 
contribution from Eq. (|15p is 5Ia oc vAB T l^hT fi,.i2- 
Consider next the reverse process, see Fig. mb), which 
yields St^^L = {tR2B{-iA)tLi) {tmtLi), and there- 
fore causes the same current contribution Sla oc {—v) ■ 
{—A)BTL,iiTji^i2. Summing up all relevant processes, 
wc find 

la oc B{tLltRl+tL2tB2){tLlAtB2+tL2{-A)tRi) (16) 

= AB [(r^ai — ri_22)rHj2 — (r^ai — TR^22)'i'L,i2] ■ 

Note that /a 7^ precisely when Eq. holds, as follows 
by explicitly computing the commutator in Eq. (|13p for 



S, QD S, S, QD S, 
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FIG. 1: (Color online) Schematic picture of transfer of a 
Cooper pair through a two-level dot. (a) Contribution to 
th^R yielding an anomalous supercurrent. (b) Reverse pro- 
cess contributing to tji-,L- Top and bottom panels represent 
initial and final states, respectively, which are connected by 
a sequence of the intermediate virtual states. Solid arrows 
indicate transitions due to tunneling (green, connecting leads 
and dot), spin-orbit (red, connecting different dot levels) and 
Zeeman (blue) coupling. For details see text. 



N = 2. The above arguments resemble the justification 
for TT-junction behavior in quantum dots with Coulomb 
blockade [l3|, and can in fact be applied to show that 
la ^ under the specified conditions even in the presence 
of interactions. Interactions give no qualitative change to 
our results since they do not alter the relevant symme- 
tries. For large Tj^^r, charging effects are smeared out in 
any case. For small Ti^/r and very strong interactions, 
the result basically equals the noninteracting one (up to 
a sign reversal of la for an occupied dot). 

Numerical analysis for N = 2. — We have numerically 
computed the full current ([5]) for a two-level dot, tak- 
ing the wave functions for a harmonic transverse and 
hard- wall longitudinal confinement as in Rcf. [lo| . De- 
noting the distance between the tunnel contacts by L, 
the Rashba coupling is kept fixed at a moderately small 
value, aL = 0.4, while the remaining parameters were 
taken both inside and outside the parameter regime 
where the restrictions needed to derive Eq. ^2]) apply. 
The resulting CPRs are shown in the main panel of Fig. [21 
where the supercurrent is plotted for increasing values of 
the Zeeman field. These plots show that a large value for 
the anomalous current la can be obtained with reason- 
able parameters. The hybridization matrices have been 
chosen to optimize the commutator in Eq. (|13p while sat- 
isfying Tl/ra2 = VTL/iJoiTTAR^- The CPR is either 
continuous or exhibits jumps associated with a change 
of the ground state, i.e., with the different occupation of 
the relevant Andreev levels. Such jumps are a common 
feature in the presence of Zeeman fields already without 
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FIG. 2: (Color online) Numerical results for the CPR of a 
two-level dot with 61 = for &^ = 0.3A, 0.7A and 0.9A (red, 
blue, brown curves, respectively, with increasing |/a|). Pa- 
rameters: aL = 0.4, Tl,11 = 2A,rL,22 = 0, Th,!! = rR,22 = 

A/2,ei = -62 = 0.65A, L = 25 nm, A = 1 meV, and 
m — 0.035me. Inset: la from numerics (red solid curve) and 
from Eq. (|12|l (black dashed curve) vs aL, with bx = 0.5A 
and same other parameters as in main panel. 



existence of such two critical currents and its potential 
for rectification behavior have been discussed in Ref. ■ 
The inset of Fig. [2] shows a comparison between the full 
numerical results for la and the analytical predictions de- 
duced from Eq. (jl2[) as a function of a. The comparison 
is done for an intermediate value of the Zeeman field, 
bx = 0.5A, and the agreement is good even for rather 
large la- 

Conclusions. — We have shown that a generic model 
for superconducting transport through a quantum dot ex- 
hibits spontaneously broken TRS, leading to an anoma- 
lous supercurrent appearing at zero phase difference be- 
tween the superconductors. The effect occurs in the pres- 
ence of spin-orbit interaction and a suitably oriented Zee- 
man field, provided that the left and right contact hy- 
bridization matrices do not commute. This implies that 
at least two dot levels must be involved. The wide avail- 
ability of mesoscopic systems holds the promise to ex- 
perimentally observe this remarkable effect in the near 
future. 



SOI Q. In the CPRs obtained by numerics from Eq. ([5]) , 
wc find that the maximal (positive) current can differ 
from the minimal (negative) current in magnitude. The 



This work was supported by the SFB TR 12 of the 
DFG, by the EU networks INSTANS and HYSWITCH, 
and by ANR-PNANO Contract MolSpintronics No. 
ANR-06-NANO-27. 



[1] A. A. Golubov, M.Yu. Kupriyanov, and E. Il'ichev, Rev. 
Mod. Phys. 76, 411 (2004). 

[2] Y.J. Doh, J.A. van Dam, A.L. Roest, E.P.A.M. Bakkers, 
L.P. Kouwenhoven, and S. De Franceschi, Science 309, 
272 (2005); J.A. van Dam, Yu.V. Nazarov, E.P.A.M. 
Bakkers, S. De Franceschi, and L.P. Kouwenhoven, Na- 
ture 442, 667 (2006); T. Sand-Jespersen et aL, Phys. 
Rev. Lett. 99, 126603 (2007). 

[3] H. Takayanagi, T. Akazaki, and J. Nitta, Phys. Rev. Lett. 
75, 3533 (1995); Th. Schapers et aL, Phys. Rev. B 67, 
014522 (2003); M. Ebel, C. Busch, U. Merkt, M. Grajcar, 
T. Plecenik, and E. Il'ichev, iMd. 71, 052506 (2005); J. 
Xiang, A. Vidan, M. Tinkham, R.M. Westervelt, and 
CM. Lieber, Nature Nanotech. 1, 208 (2006); S. TireUi, 
A.M. Savin, C.P. Garcia, J. P. Pekola, F. Beltram, and F. 
Giazotto, Phys. Rev. Lett. 101, 077004 (2008). 

[4] For recent work and references to previous nanotube 
experiments, see A. Eichler et aL, Phys. Rev. B 79, 
161407(R) (2009). In nanotubes, the effects described 
in this paper are expected to be very small because the 
SOI is weak and the tunnel contacts often preserve valley 
pseudospin. 

[5] R. Winkler, Spin-Orbit Coupling Effects in Two- 
Dimensional Electron and Hole Systems (Springer, New 
York, 2003). 

[6] J. Nitta, T. Akazaki, H. Takayanagi, and T. Enoki, Phys. 

Rev. Lett. 78, 1335 (1997); G. Engels, J. Lange, T. 

Schapers, and H. Liith, Phys. Rev. B 55, R1958 (1997). 
[7] E.V. Bezuglyi, A.S. Rozhavsky, I.D. Vaguer, and P. 



Wyder, Phys. Rev. B 66, 052508 (2002). 

[8] I.V. Krive, L.Y. Gorelik, R.I. Shekhter, and M. Jonson, 
Fiz. Nizk. Temp. 30, 535 (2004) [Low Temp. Phys. 30, 
398 (2004)]; I.V. Krive, A.M. Kadigrobov, R.I. Shekhter, 
and M. Jonson, Phys. Rev. B 71, 214516 (2005). 

[9] O.V. Dimitrova and M.V. Feigel'man, JETP 102, 652 
(2006). 

[10] L. DeirAnna, A. Zazunov, R. Egger, and T. Martin, 
Phys. Rev. B 75, 085305 (2007). 

[11] B. Beri, J.H. Bardarson, and C.W.J. Beenakker, Phys. 
Rev. B 77, 045311 (2008). 

[12] A. A. Reynoso, G. Usaj, C.A. Balseiro, D. Feinberg, and 
M. Avignon, Phys. Rev. Lett. 101, 107001 (2008). 

[13] A. Buzdin, Phys. Rev. Lett. 101, 107005 (2008). 

[14] V.B. Geshkenlsein and A.I. Larkin, Pis'ma Zh. Eksp. 
Teor. Fiz. 43, 306 (1986) [JETP Lett. 43, 395 (1986)]; 
S. Yip, Phys. Rev. B 52, 3087 (1995); Y. Tanaka and S. 
Kashiwaya, ibid. 56, 892 (1997); Y. Asano, Y. Tanaka, 
M. Sigrist, and S. Kashiwaya, ibid. 71, 214501 (2005). 

[15] The prediction /„ / in Ref. [H is an artefact of an 
approximation violating a gauge invariance of the sys- 
tem. This problem can be traced to an inconsistent use 
of the quasiclassical approximation. For a correct treat- 
ment, one should proceed along the lines of Ref. 0]. 

[16] The difference to Eq. ([2| arises because a gauge transfor- 
mation is necessary to obtain the exact solution [lo| . 

[17] B.I. Spivak and S.A. Kivelson, Phys. Rev. B 43, 3740 
(1991). 



